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Abstract 

We show that many two dimensional domain patterns observed in Monte 
Carlo simulations can be obtained from the many soliton solutions of the 
imaginary time Sine Gordon equation. This opens the door to analytic phys- 
ical understanding of the micromagnetics in ultra-thin films. 
There has been much experimental interest recently in the magnetism of ultra-thin films. 

partly motivated by the possible integration of the semi-conductor microelectronics 
technology with magnetic elements [0 and possible device applications with the giant mag- 
netoresistive (GMR) effect. From a fundamental physics viewpoint, these systems present 
opportunities for studying new phenomena that are beginning to be uncovered. The inter- 
action energy between the spins at positions R, R' is 

# = 0.5 ]T V^R-RO^R^R') (1) 

ij=xyz,KR' 

where V = Vd + V e + V a is the sum of the dipolar energy Vdjj(R) = gVjVj(l/|R|); the 
exchange energy V e = — J5(R = R' + d)^-; and the crystalline anisotropy energy V a . Here d 
denotes the nearest neighbours, g and J are coupling constants. The form of the anisotropy 
energy depends on the material of interest. It can be uniaxial (e.g. V a = —KJ2iSf x ) or 
four-fold symmetric (e.g. V a = —KJ2i[Si X ~ ^iy] 2 /^)y with the easy or hard axis aligned 
along specific directions. The dipolar interaction often lead to the formation of domains. 
The pattern of the domains have recently received considerable interest under the context 

1 



of the "self-assembled" systems where electric dipoles lead to the formation of domains 
in Langmuir films. Whereas the electric dipoles are always perpendicular to the film plane 
in that case, the magnetic dipoles can be parallel or perpendicular to the plane For 
discussions in this paper, we restrict our attention to those cases so that the spins lie in the 
plane of the film, the case of experimental interest in sensor type applications. 

The domain pattern depends on the shape of the sample, which is especially important 
for small structures. The physics of the pattern of domains in small magnetic structures is 
the subject of the present paper. We have been studying the physics of spin reversals of 
different small structures such as monolayer films with perpendicular || and 4- fold in 
plane || anisotropy, nanowires and particles fLOf , coupled films |L1| and the shape of the 



nucleus |TI|. This paper reports our findings that much of the domain patterns observed 



in the numerical simulations can be reproduced as the analytic many soliton solutions 
of the imaginary time Sine-Gordon equation. This is illustrated by two examples in Fig 1 
and 2 where we show the simulation and analytic results side by side. These analytic results 
have the potential for greatly improving understanding quantitatively the domain structure 
and the switching process in small structures. Thus analytic calculations can be performed 
to predict trends as the system parameters are changed. These analytic results can be used 
as a starting point of a simulation, considerably shortening the simulation time; sometimes 
the simulations become entirely unnecessary. We now explain our results in detail. 

Mathematically in the continuum approximation, the dipolar energy Ed ~ 
| / c?Rc/R / S'j(R)S'j(R / ) VjVj(l/|R — R'|) can be written in terms of the magnetic charges 
V • S after two integration by parts and neglecting the boundary terms as Ej « | / c/Rc/R'V • 
S(R)V • S(R')(1/|R — R'|). Thus the dipolar energy is reduced if the "magnetic charges" 
are as small as possible. This is usually achieved when lines of dipoles form closed loops. 
The orientation of the spin is determined by its angle <fi. For example, when the azimuth 
angle 9 can be described as a vortex with = 9 — n/2 the dipolar energy is minimized. 
When this type of global constraint is satisfied, the domain structure is usually determined 
by minimizing the exchange and the anisotropy energy; we obtain the equation: 



O.5ATsin40/J = 0, 



(2) 



Here J ~ zJ/A is the effective exchange, z is the number of nearest neighbours. It comes 
from converting the discrete model to the continuum approximation. The exactly soluable 
sine-Gordon equation — d 2 )<f) — O.5i^sin40/J = is formally the same as the above 
equation (0) if we transform the y coordinate into the imaginary time it. In this way, we 



can generate a 90 degree domain wall "soliton" solution as <fi = tan 1 exp[— J2K/Jx] where 
the angle <ft changes by 90 degree as the wall is traversed and x changes sign. This solution 



is one dimensional and is well known p"3 |. 

Many soliton solutions are known but have never been exploited in the understanding 
of domain structures. A general two soliton solution of the sine-Gordon equation has the 
form ( 0): 
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Here x' = xy2K/J,t' = ty2K/J. This solution has 4 arbitrary constants Ui, U2, x\, x%. 
Using the transformation u\ = iv,u,2 = —iv,t' = i(y' — y'x)-,x' x = x' 2 and choosing x[ = 
(ln(l/t>) —i7r/2)/j + x' , y[ = —in/ (2vy) +y' , we obtain the two soliton solution in the form: 



= tan" 1 [sinhfrv J2K/J(y - y ))/(-v sinh( 7 J2K/J(x - x )))] 



(4) 



where v is a parameter, 7 = 1/y/l + v 2 . This describes a closure domain. An example is 
shown in Fig. IB for a triangular lattice of 3600 spins for K = 0.2 and J = 2. A closure 
domain can be viewed as the space-time trajectory of two solitons coming together and 
eventually moving apart. The parameter v describes the orientation of the domain wall. 
To fit into a sample of aspect ratio r, one expects v — r, as we have verified directly by 
numerical calculation. For a triangular lattice, the center of the defect, (xo,y ) for the 
lowest energy configuration sits in the middle of the triangle. This type of domain wall is 
often observed in simulations in systems in zero external magnetic field. A typical finite 
temperature simulation result [jllj is also shown in Fig. 1A for the same value of J and 



K and g = 1, obtained from cooling a high temperature configuration that starts off with 
all spins aligned in the x direction. To study the possible effect of the dipolar interaction 
and the accuracy of the analytic formula, we have numerically minimized the total energy 
of the system starting from the configuration given by the analytic formula and using a 
quasi-Newton algoraithm for a system with 400 spins. We have explored different values 
of g less than 1 and find that the mean square difference between the initial and final 
azimuthal angles is less than 0.1 radian, out of a possible range of tt. Thus the accuracy 
is 3%; the analytic formula is indeed a good approximation. With this analytic formula, 
it is much easier to investigate the physical properties of closure domains quantitatively. 
For example, we have investigated the size dependence of the energy difference between the 
closure domain and that with uniform magnetization along the x direction. The difference 
in energy divided by the effective coupling constants (g for the dipolar energy and yJK 
for the sum of the exchange and the anisotropy energy) is shown in Fig. 3 as a function 
of the sample size. For a rectangular sample of a triangular lattice with an aspect ratio of 
0.866 and x dimension L\, the dipolar energy difference AE p can be fitted by the formula 
g(109.5 — 10.54Li) with an error of less than 4% whereas the sum of the anisotropy and 
exchange energy AE W can be fitted by the formula \J JK(28.76 + 2.641/!) with an error 
of less than 0.3%. The closure domain is lower in energy than the uniformly magnetized 
state when the sum of these two energies become negative. For a film of thickness d, we 
expect that approximately g = god 2 , J = J d, K = K d where the subscript refers to 
the bare coupling per spin. Thus the closure domain is lower in energy for sample sizes 
Li > (109.5 + 2.876 v / Jo^o/#o<)/(10.54 - 2My/J^/g d). This can only happen if the 
denominator is positive; ie d > d c = 0.25^ JqKq/ g^. As an example, consider bcc Fe where 
go ~ 0.254i^, K « 0.038/T, Jo ~ 500K. Thus the closure domain is lower in energy for 
thicknesses d > 4.3 layers. 

The solutions (f|) is, strictly speaking, applicable to infinite samples. The consideration of 
the domain patterns in small structures require the imposition of finite boundary conditions. 
The solutions of Eq. (fj) which satisfy these boundary conditions can be obtained starting 

4 



from the ansatz suggested by Lamb [15| for the solution of the sine-Gordon equation. We 
seek solutions of Eq. (|2|) having the form <p(x,y) = tan~ 1 [/ (x')g(y')], where / and g are, 
in general, Jacobian elliptic functions defined by [[Uj (f) 2 = a/ 4 + j3f 2 — 7 and (g') 2 = 



— 7(7 4 — l)g 2 + a with a, (3, and 7 arbitrary constants, x' = y2K/ Jx, y' = y2K/ Jy. As 
an example, we consider configurations corresponding to edge domains with the boundary 
conditions that the spins point up (down) on the left (right) edge and horizontally on both 
the top and the bottom edge. 

cn^l + k 2 y', ki g ) 



tan 1 



Atn(Qx', X 



(5) 



dn(tyi + k 2 y\k lg 

where k 2 = [A 2 {l 2 (l - A 2 )]/[{l 2 (l - A 2 ) 2 - 1], k\ g = A 2 tl 2 (l - A 2 )/(tt 2 (l - A 2 ) - 1), 
X 2 = [A 2 + n 2 (l - A 2 ) 2 ]/[n 2 (l - A 2 )} and v 2 = [fi 2 (l - A 2 ) 2 - 1]/[1 - A 2 }. The parameters 
A and Q can be determined by requiring that the component of S normal to the surface 
boundary be zero so that the dipolar energy is minimized. 

Figure 2B shows the edge domain pattern obtained by using Eqs. ([5]) for a triangular 



lattice 3600 spins for J = 2 and K = 0.2. In Figure 2A we show the Monte Carlo result [TT 
for a bilayer system for a triangular lattice of 3600 spins for the same value of J and K and 
g=l. Similar domain patterns are also seen in the zero field remanent state for a system 
with a single layer. || 

To study the possible effect of the dipolar interaction and the accuracy of the analytic 
formula, we have numerically minimized the total energy of the system starting from the 
configuration given by the analytic formula and using a quasi-Newton algoraithm for a 
system with 400 spins. When the dipolar interaction is too small, our algorithm recovers 
the minimum energy state of uniform magnetization. As long as the dipolar interaction is big 
enough the minimum energy configuration from our algoraithm is essentially independent 
of the strength of the dipolar interaction. We obtain a state that resembles our analytic 
results. The mean square difference between the initial and final azimuthal angles is less 
than 10%. The analytic formula is indeed a good approximation, even though it is not as 
good as that for the closure domains. With this analytic formula, we have investigated the 



size dependence of the energy difference between the edge domain and that with uniform 
magnetization along the x direction. The results are shown in Fig. 3. For a rectangular 
sample of a triangular lattice with an aspect ratio of 0.866 and x dimension Li, the dipolar 
energy difference AE P can be fitted by the formula g(52.87 — 3.97Li) with an error of less 
than 5% whereas the sum of the anisotropy and exchange energy AE W can be fitted by the 
formula y/ JK(10A6 + 1.9Li) with an error of less than 1%. The edge domain is thus of lower 
energy when the sum of these two energies become negative. As expected, when compared 
with the closure domains, the dipolar energy gained is less while the cost in the anisotropy 
and exchange energy is also smaller. For a film of thickness d, the edge domain is lower in 
energy for sample sizes L\ > (52.87+ 10. 46 \J JqKq/ god) / (3.97 — 1.9\/ JoK / g d). This can 
only happen if the denominator is positive; ie d > d ec = 2\J JoK /g . For bcc Fe, d ec = 8.2 
layers. 

In this paper we have discussed two examples of analytic solutions for domain patterns. 
Many possibilities remain to be explored. For example, consider 

(p = ta,n- 1 [cos(- i v^2K/J(y - yQ)) / {v swh{<i \] 2K / J {x - x ))]. (6) 

where 7 = 1/yl — v 2 . This solution can be considered the analytic continuation of the 
solution (§) with an imaginary v. When v is small, this solution describes two 90 degree 
domain walls separated by a distance 2 \n(2/v)\fj/2K/j). As v is increased from zero, two 
separated 90 degree domain wall merge to become a 180 degree domain wall with vortices in 
between. This type of solutions is not the lowest energy configuration in zero magnetic field 
but occurs as a rate limiting step in spin reversal processes at a finite magnetic field. Our 
solution provides for configurations that are local extrema of the exchange and anisotropy 
energy. The ordinary 180 degree domain wall in zero field, which conisists of two 90 degree 
domain walls, is not a local extrema of the exchange and anisotropy energy. It is only 
stabilized by the magnetoelastic or dipolar energy. |J 

In summary, we have provided examples of how the many soliton solutions can be used 
to understand the domain structures in ultra-thin films. This opens the door to analytic 
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quantatitive understanding of the micromagnetics in these systems. 
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FIGURES 

FIG. 1. Closure domain configuration for a rectangle from a 2 soliton solution. The analytic 
results is in B. The finite temperature Monte Carlo results observed in ref. 11 is shown in Fig. A. 

FIG. 2. Edge domain configuration for a rectangle from a 2 soliton solution. The analytic 
results is in B. The finite temperature Monte Carlo results observed in ref. 11 is shown in Fig. A. 

FIG. 3. The energy difference between the domain configuration and that with uniform magne- 
tization as a function of the linear dimension of the sample. These energy differences are normalized 
by the coupling constants, as is described in the text. 
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